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Abstract. We develop a microscopic theory of electron transport in superlattices 
within the Wannier-Stark approach by including the interaction associated with Zener 
tuneling among the energy levels pertaining to adjacent quantum wells. By using 
a Monte Carlo technique we have simulated the hopping motion associated with 
absorption and emission of polar optical phonons and determined the main transport 
parameters for the case of a GaAs/GaAlAs structure at room temperature. The 
interaction among the levels is found to be responsible for a systematic increase of 
the level energy with respect to the bottom of the quantum well at electric fields above 
about 20 kV/cm. When compared with the non-interacting case, at the highest fields 
the average carrier energy evidences a consistent increase which leads to a significant 
softening of the negative differential value of both the drift velocity and diffusivity. 
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1. Introduction 



In recent years, superlattices (SLs) received a relevant scientific and technological 
interest, owing to their nonlinear electrical and optical properties. In particular, at 
increasing electric field strengths SLs exhibit a nonOhmic behavior with the presence 
of a negative differential conductivity (NDC) region [1] above a threshold field. In 
the nonOhmic regime, charge transport can be described in terms of hopping between 
Wannier-Stark (WS) states E], thus SLs are good systems for studying hopping 
transport in semiconductor structures. Despite of the above interest, the knowledge of 
the transport parameters, like drift velocity and diffusion coefficients for a given SL is 
still mostly an open problem. 

From the analytical side, the problem was investigated with different approaches 
in Refs. [[U El EJ [7J IB]]. The coupling between energy levels and the electric field (the 
so called Zener effect [9]) was found responsible for the onset of a resonance in the 
drift velocity [HE] and of an antiresonance in the diffusion|6] at fields where both these 
quantities are decreasing functions of the applied electric field. 

From the simulative side, recent Monte Carlo investigations [TUl ITT] have analyzed 
the drift velocity [10] and both the drift velocity and diffusion [TT] but neglecting the 
coupling between the field and the energy levels. Furthermore, in Ref. [ |10| ] only general 
trends are reported so that a quantitive evaluation of the drift velocity is absent. 

Here we develop a microscopic theory of electron transport in SLs within the 
Wannier-Stark approach, which includes the coupling among the energy levels due to the 
presence of a strong electric field so as to cause Zener tunneling. The theory is applied 
to an ideally infinite GaAs/GaAlAs superlattice at 300 K. The electric fields considered 
are in the range 5 -j- 200 kV/cm where the Wannier Stark approach is well justified 
[T2| [13], for which experiments are available |14j . A comparative investigation of the 
effect of the interaction among energy levels with respect to the case of independent 
energy levels is then carried out. 

The content of the paper is organized as follows. In Sec. 2 we define the physical 
system under investigation and present the theoretical approach to be used within a 
Monte Carlo simulator. In Sec. 3 the main results concerning the drift velocity and the 
diffusion coefficient are presented. Major conclusions are drawn in Sec. 4. 



2. Theoretical approach 

2.1. Theory of WS states 

We consider a SL made by alternated layers of different semiconductors (A and B) along 
the ^-direction. The single electron Schrodinger equation in the presence of an external 
electric field along the ^-direction F = (0, 0, F) is: 

: 2 
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2m* 



V sl (t) - eFz 



(1) 
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Figure 1. Eigenvalues of the first and second energy levels of a quantum well 
pertaining to an ideal infinite SL. Symbols refer to the case of interacting energy 
levels, dashed lines to the case of independent energy levels. 



where p is the momentum operator, m* is the effective mass assumed to be the same in 
the well and in the barrier, Vs%(r) is the potential which determines the barrier height 
of the SL, and ip(r) = (p(z)g(x,y) with g the Bloch function. The SL potential is taken 
of the Kronig-Penney form: 

. v / V if < z < a 
Vsl{z) = \0 iia<z<a + b (2) 

with a and b the length of the well and of the barrier, respectively, and d = a + b 
the SL period. Accordingly, using the Dirac notation the Schrodinger equation in the 
^-direction is: 

\H -eFz]\tp)=e\<p). (3) 

The solution of Eq. (j3J) can be expanded over the SL Bloch states \n,k z ), with n the 
band index, which are eigenstates of the Hamiltonian Hq of the SL in the absence of 
the electric field 

l?>=£ I dk'Mn'l^lKn') (4) 

n! z 

and substituted into the above Schrodinger equation; multiplying by (k z n\ we develop 
the calculation and obtain 

[E n (k z ) - e] (k z n\ip) 

eF^2 I dk' z (k z n\z\k' z n')(k' z n'\<p) = . (5) 
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Figure 2. Normalized square value of the Wannier Stark wavefunction at different 
electric field strengths of, respectively from right to left, 5, 30 and 200 kV/cm. The 
first raw of figures refers to the first level in the absence of interaction; the second raw 
of figures refers to the first level in the presence of interaction; the third raw of figures 
refers to the second level in the presence of interaction. 



Here, the value (k z n\z\k' z n') is calculated explicitly and, after some algebra, it results to 
be 

(k z n\z\k' z n f ) = ^r-^ 7 5 nn >5{k z - k' z ) + X nn ,5{k z - k' z ) (6) 



X nn , = ^ / &zu* nkz (z)-^-u n , K (z) . (7) 



with 

X , _ : 

L J n/CzV J dk' z 

where u n k z is the periodic part of the Bloch function. The above matrix elements include 
the interaction among levels pertaining to different quantum wells caused by the electric 
field and they are basically related to Zener interband tunneling. We recall, that under 
the conditions of moderate electric fields and weak coupling between adjacent wells it is 
justified to consider a one-band model for the calculation, thus neglecting off-diagonal 
matrix elements described by eq. (JTj). This approximation leads to a simple solution of 
eq. ([5]) and was already applied in Ref. [[10J [11]] where the WS eigenstates have been 
obtained analytically [TT] . 
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2.2. Numerical method 



In the present paper we point to fully take into account the off-diagonal terms (|7j) in the 
calculation of the WS eigenstates and further of the SL transport properties. However, 
in this case it is not possible to get a simple analytical solution of the system. Thus, to 
accomplish this purpose, we have decided to solve numerically the Schrodinger equation 
along the z direction ([3]). With this approach, no interaction is dropped during the 
calcualtion and the interaction decribed by (171) is automatically included in the solution. 
Looking for a solution of the Schrodinger equation at the finite differences, we consider 
a number M of cells (typically 30 cells in the backward direction and 30 in the forward 
direction for a total number of M = 61) representative of an ideal infinite SL. (Indeed, 
electron wavefunctions are found to extend up to a maximum of 20 cells). Then, the 
Schrodinger equation is discretized in a number of N points (typically N= 100) equally 
spaced inside a single cell, M x N being the total number of points used to obtain 
the numerical solution within an estimated convergence of a few percent. Accordingly, 
z — ► £A, <f(z) — > cfi, with A = 4, and defining 
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the discretized Schrodinger equation of the system takes the compact form: 

Aip = sf (8) 

with if = (<f i, if 2 , ■ ■ -,<Pnm)- 

Since the matrix A is real and symmetric the eigenvalues are real, distinguishable, 
and in number equal to the order of the matrix M x N. The eigenvectors f % are 
the Wannier-Stark wavefunctions. In the present case, the (fj are two-dimensional 
vectors with the index % representing the corresponding energy level and the index j 
the pertaining grid point. The eigenvalue e l and the eigenvector f l , being labelled with 
the same index i, are in full correspondence. 

The eigenvalues are then classified according to the pertaining quantum well and 
energy level by using the properties of the Wannier-Stark eigenvectors (for a review of the 
WS functions properties see for example |15j). In particular we know that WS functions 
should be localized and satisfy the translational symmetry f n+m (z) = f n (z — md) 
for any pairs of quantum wells labelled by (n,m). Then, as validity test of the 
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numerical algorithm, the energy difference between adjacent wells is calculated and 
their ladder distribution verified. We have also verified that in the central region of the 
representative device, where boundary effects are negligible, the translational property 
of the eigenvectors belonging to the same energy level is accomplished. In this way, for 
each value of the electric field we select the wavefunction representative of the infinite 
SL. 



2.3. Electron-phonon interaction 



The calculation of the polar-optical phonon scattering probability per unit time in the 
WS representation follows the standard Fermi Golden Rule. In brief, after decoupling 
the phonon bath, the point-to-point transition probability per unit time is 



P(unk\\, v'n'\s!w 



h ^ 

it, q 



c 2 (q) 



x 



(v'n'\e TiqzZ \isn) 



5, 



k ii- k ([=Fqi 



5 



e(nki|) — e(n'kn) =f hu op + veFd 



(9) 



where n q is the mean value of the phonon occupation number, q = (qn, q z ) is the phonon 
wave vector and k = (k\\,k z ). For c 2 (q) we take the standard Frolich squared matrix 
element 1161 : 

(10) 



c 2 (q) 



he 2 uj op 



1 



1 



q 2 ,) 2 



2Ve \e + e_J (q 2 

where u op is the optical phonon frequency, e_, e + are the static and high frequency 
relative dielectric constant, q D is the screening wave vector and V is the crystal volume. 

Being interested in the total scattering probability per unit time from an initial state 
(imk[|) to any possible final state kn, we integrate over the final momentum variable and 
obtain the integrated transition probability per unit time: 



d6 



m 



n 

n + 1 



' (2tt) 5 ^ 3 
dg 2 c 2 (q) (v'n'\e Tig * z \vn) 



x 



Having integrated the 5 function in Eq. ([9]) implies having accounted for momentum 
conservation that leads to the constraint q = (k' x — k x , k' — k y , q z ). This last Eq. ffTTl) 
is not further integrable with standard method because of the not simple dependence 
upon q. Within the Monte Carlo technique this step is accomplished using the rejection 
technique [T7] . 

Equation (lllj) describes the probability that an electron with initial state (Vnkii) 
is scattered to a final state in well u', in band n', with any ky. However, in our 
approach electrons are allowed to jump a maximum distance of 6 wells forward and 
6 wells backward, since the hopping transition element is negligible for longer distances. 
Indeed, we have verified that, for the transport parameters investigated here, the results 
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are not affected by a further increase of the maximum hopping distance. Accordingly, 
the associated transition probabilities per unit time for a carrier jump from an initial 
well to any of the 12 final wells due to absorption and emission of polar optical phonons 
is evaluated and used inside a Monte Carlo procedure, as detailed in Ref. [[TT]]. For 
sake of simplicity, other scattering mechanisms are neglected and we have considered 
at most only the first two energy levels in each well. The time evolution of the carrier 
ensemble is then evaluated for a period sufficiently long to extract the drift velocity and 
the longitudinal diffusion coefficient from the corresponding slopes of the average value 
of the position and of its variance along the field direction [TT] . 

3. Results and discussion 

The theory is applied to the case of a GaAs/GaAlAs SL at 300 K, being that the most 
studied SL. As typical parameters we take a = 8.5 nm, b = 1.5 nm and a barrier 
height of 0.4 eV. For the electron effective mass we take everywhere the value of GaAs 
m* = 0.067 mo, with itlq the free electron mass, for the polar optical phonon energy 
we take e op = 35 meV, and for the polar optical coupling we take — — — = 0.013. 
The electric fields considered are in the range 5 4- 200 kV/cm where the Wannier Stark 
approach is well justified [T3] . 

An ensemble Monte Carlo simulation has been set-up in order to simulate electron 
transport transport in this superlattice, as detailed in [[TT]] and the estimated physical 
quantities are calculated as follows 

1 N 

D. = (13) 
2 6t K ' 

1 N 

<*> (14) 
i=i 

for drift velocity, diffusion coefficient end mean energy respectively, where Zi and q are 
the position and the energy of the i-th particle. 

The results of the investigation are reported in Figs. 1 to 5. Here, the case of 
interacting levels (through Zener tunneling) is compared to the case of independent 
levels at increasing values of the electric field. 

Figure 1 reports the energy of the first and second level inside the quantum well 
(measured from the bottom of the well) as function of the applied field. In the absence of 
interaction the energy of the first and second levels are found to be of 43 and 175 meV, 
respectively. The presence of the interaction is responsible for a significant increase, up 
to about a factor two, of the energy levels at fields starting from about 10 kV/cm. Such 
an increase is associated with the broadening of energy levels [IEJ [T9j . 

Figure 2 reports the square modulus of the Wannier Stark wavefunctions associated 
with the first and second energy level at increasing electric fields (from left to right). The 
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Figure 3. Average carrier energy as function of the applied electric field. Open 
circles, full triangles, and open squares refer to the cases when only the first level in 
the absence of interaction, only the first level in the presence of interaction, both the 
first and second level in the presence of interaction are considered. Curves are guides 
to eyes. 



set in the first (second) raw of the figure refers to independent (interacting) first energy 
levels. The tendency to localization exhibited by both cases at moderately increasing 
fields is found to be broken at the highest fields by the presence of the coupling between 
level and field, as evidenced by the long tail in the direction of the field exhibited by the 
wavefunction in the second raw. This tail is originated by the coupling of different bands 
of the superlattice induced by the electric field, as described in Eq. ((?]), and enables us 
to evaluate the effect of carrier tunneling from a well to those in the forward direction 
(Zener effect [9]). The agreement between analytical and numerical results at moderate 
electric fields is taken as a validation test of the numerical algorithm used here. The set 
in the third raw of Fig. 2 refers to the case of interacting second energy levels. Here, 
at low fields the derealization is stronger than in the case of the first level, as expected 
for a higher energy level. At increasing fields we observe an increase of localization and 
the appearance of the long tail in the direction of the field. 

Figure 3 reports the carrier average energy obtained from Monte Carlo simulations 
at increasing electric fields. Above 50 kV/cm, the presence of interacting energy levels 
is found to be responsible for a substantial increase of the average energy (up to about 
a factor of 10) associated with the increased spatial mobility of the carriers due to the 
Zener effect. The inclusion of the second energy level is found to emphasize the trend 
of energy increase.. 
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Figure 4. Average carrier velocity as function of applied electric field. Full triangles, 
open circles and open squares refer to the cases when only the first level in the absence 
of interaction, only the first level in the presence of interaction, both the first and 
second level in the presence of interaction are considered. Curves are guides to eyes. 



Figures 4 and 5 report the average velocity and the longitudinal diffusion coefficient 
at increasing fields. Apart from a numerical off-set within a factor of two, in the 
region up to about 50 kV/cm, both quantities evidence the same systematic decrease 
with increasing fields. Such a decrease is associated with the localization trend of the 
wavefunction [20] and it is called negative differential mobility (NDM). Above 50 kV/cm, 
the simulations show a substantial increase (up to about a factor of 10) of both drift and 
diffusion with respect to the independent case, in concomitance with that of the average 
energy This softening of the negative diferential behaviour of both drift and diffusion is 
associated with the onset of the Zener effect. The inclusion of the second energy level is 
found to emphasize the trend to increase of both the kinetic coefficients. We remark that 
simulations evidence an appreciable kink of both drift and diffusion at a field around 
40 kV/ cm, which is associated with the resonance occurring when the potential energy 
due to the voltage drop on a period equals the optical phonon energy [U ITT] . 

For a microscopic interpretation of the above results, we refer to the shape of the 
total hopping rate from an initial well v to other final wells v' for a forward or backward 
hop in the simple case of a single level, P uu > . The maximum value of these rates, which 
are used in Monte Carlo simulations to generate a hop between states, are reported in 
Figs. 6 to 8 at different field strengths for the case of independent and interacting first 
energy levels. 

Figure 6 shows that at low fields the independent and interacting models give 
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Figure 5. Carrier longitudinal diffusion coefficient as function of the applied electric 
field. Full triangles, open circles and open squares refer to the cases when only the first 
level in the absence of interaction, only the first level in the presence of interaction, 
both the first and second level in the presence of interaction are considered. Curves 
are guides to eyes. 



practically the same results (we notice that the forward value is always greater than the 
corresponding bacward value, thus ensuring a net motion of the carrier ensemble along 
the field direction). Hopping rates at increasing distances decreases exponentially as 
expected by the localized feature of the electron wavefunction. Figures 7 and 8 show 
that the presence of the interaction leads to a systematic increase of the rates expecially 
of those to far neighbour wells. Accordingly, the drift and diffusion in the presence of 
interacting levels increases with respect to the case of independent levels. 

The general behaviour of the drift velocity is in good qualitative agreement with 
the analytical results of Ref. [[3]] in the common range of fields. Furthermore, recent 
experiments [H] on the drift velocity in GaAs/GaAlAs SLs evidence a saturation 
behavior of the drift velocity for electric fields in the range 20 -j- 50 kV/cm, which is in 
qualitative agreement with present results. In this respect, we notice that for a proper 
experimental validation of the theory more quantitative experiments are necessary. In 
this context, we remark that present results on the diffusion coefficient do not evidence 
the Zener antiresonance predicted by analytical calculations [6] . The lack of a sufficient 
resolution of the results due to the numerical uncertainty, estimated to be within a 
factor 2 at worst, and/or the presence of a substantial carrier heating at the highest 
fields, could be the reasons for the missed evidence of the antiresonance in the diffusion 
coefficient. 
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Figure 6. Maximum hopping rate from a given inital well v = to a final well in the 
forward (positive) and bacward (negative) direction. Calculations are performed for 
the case of a single level and E = 13 kV/cm. Full circles refer to independent levels, 
full triangles to interacting levels. Curves are guides to eyes. 
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4. Conclusions 

We have investigated transport properties of an ideal infinite superlattice in the presence 
of interacting energy levels associated with the presence of an applied electric field. To 
this purpose the one dimensional Schrodinger equation in the direction of the applied 
field has been numerically solved and an ensemble Monte Carlo technique used to 
simulate the carrier motion in the three-dimensional states of a perfect superlattice. 
Theory is applied to a GaAs/GaAlAs structure with a period of 10 nm at room 
temperature. The main results of the present work can be summarized as follows, (i) the 
interaction among the levels is found to be responsible for a systematic increase of the 
energy levels with respect to the bottom of the quantum well, (ii) the same interaction 
is found to be responsible for the onset of a new derealization of the Wannier-Stark 
states at the highest fields associated with the tunneling of electrons between adjacent 
quantum wells, (iii) with respect to the non-interacting model, for fields above about 
20 kV/cm the average carrier energy as well as the drift velocity and the longitudinal 
diffusion coefficient evidence a consistent increase (up to about a factor of 10). Such 
an increase leads to a significant softening of the negative differential value of both the 
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Figure 7. Maximum hopping rate from a given inital well v = to a final well in the 
forward (positive) and bacward (negative) direction. Calculations are performed for 
the case of a single level and E = 50 kV/cm. Full circles refer to independent levels, 
full triangles to interacting levels. Curves are guides to eyes. 




drift velocity and diffusivity at fields where Bloch oscillation regime should be active. 
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